Nonequilibrium phase transitions in active contractile polar filaments 
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We study the patterning and fluctuations of a collection of active contractile polar filaments 
on a two dimensional substrate, using a continuum description in the presence of athermal noise, 
parametrized by an active temperature Ta ■ The steady states generically consist of arrays of inward 
pointing asters and show a continuous transition from a moving aster street to a stationary aster 
lattice. In contrast to its equilibrium counterpart, this active crystal shows true long range order 
at low Ta- On increasing Ta, the asters remodel with a distribution of lifetimes; concomitantly we 
find novel phase transitions characterized by polar and bond-orientational order. 



A remarkable feature of living cellular systems is that 
the same evolutionarily conserved ingredients - filaments 
(actin), motors (myosin) and their regulators, in the pres- 
ence of ATP - can exhibit a variety of phenotypes de- 
pending upon the cellular context pQ. Recently, there 
have been serious attempts to quantitatively understand 
this spectrum of behaviours using in-vitro reconstituted 
systems - In this paper, we provide a detailed anal- 
ysis of the nonequilibrium steady states and phase tran- 
sitions of actively driven collections of filaments-motors 
in two dimensions using the framework of active hydro- 
dynamics or active gels [4,5, which could in principle be 
used to compare with such experiments. 

Our study differs from earlier studies of active hydro- 
dynamics in two important aspects : (i) We provide 
an exhaustive treatment of the steady states and their 
transitions using dynamical equations written in terms 
of both the concentration and polarization of the active 
filaments. Realizing that there are many different micro- 
scopic processes, such as acto- myosin contractility and 
treadmilling, that may simultaneously engage with the 
cortical actin, has allowed us to probe parameter regimes 
that have not been explored before, (ii) We study the ef- 
fects of spatiotemporal active noise that are inevitably 
present in cellular systems. This not only affects the dy- 
namics in the steady state, but also induces novel phase 
transitions that are characterized by a variety of order 
parameters. Our method of analysis combines both an- 
alytical (linear stability analysis) and detailed numerical 
solutions and complements [5]. Our main results: (a) 
We find a variety of steady state configurations that in- 
clude domain walls, boojums, inward-pointing asters and 
spirals; vortices are generically unstable. The phase di- 
agram, which includes a transition from a moving aster 
street phase to a stationary aster lattice phase across a 
critical line (Fig. 1), is robust and is our first main re- 
sult, (b) The 2d active aster lattice phase is stable even 
in the presence of noise; it shows true long-range order 
(LRO), in striking contrast to its equilibrium counterpart 
(Fig. 2). (c) Beyond a critical active temperature, this 
transforms to an aster lattice with quasi long range or- 
der (QLRO, Fig. 2) during which the asters remodel with 
a power law distribution of lifetimes. This phase exhibits 
strong bond tetratic order (Fig. 3). (d) On increasing 




Figure 1: (a) Phase diagram in K2/K1 versus C/ v o at Ta = 0, 
where D = 2.5 and a = f3 = 100 in dimensionless units 
(see text), shows three phases - (i) Boojum stripes (ii) Aster 
street and (iii) Aster lattice. Representative configurations of 
c (shading) and n (arrows), when K\ = 2.5, K2 = 0, Wq = 1, 
and £ = 0.1 (i), 5 (ii) and 80 (iii) are shown alongside, (i) and 
(ii) are moving phases, while (iii) is a stationary phase, (b) 
The net polarization (en) goes to zero continuously, as the 
contractility parameter C/«o changes across the moving to 
stationary phases, with a power law behaviour (en) ~ [£ — 
Cc| 7 , where £ c = 40 and 7 ~ 1.4 (Inset shows a logarithmic 
plot (base e), lines show fit). 



the active temperature further, we observe a discontinu- 
ous transition to a bond nematic phase (Fig. 3) with an 
exponential aster lifetime distribution. Significantly, this 
bond nematic phase is polar, the filaments have a net ori- 
entation and exhibit power-law orientational correlations 
in this phase (Fig. 4), with exponents in the Toner- Tu 
universality class 

We describe the acto-myosin system as a collection of 
permanent force-dipoles, which being polar, will induce a 
net drift of the i th filament with respect to the medium, 
represented by a polarization vector n^. The hydrody- 
namic fields are the local concentration c(r, t), the polar- 
ization density c(r, t)n(r, t), and the hydrodynamic ve- 
locity v(r,t) [2H9]. 

We will assume that momentum is lost by local friction 
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at the "substrate", thus Tv = —V ■ a, where T is the 
friction coefficient and a oc cnn is the active stress (TO] 
due to the force-dipoles. We can use this to eliminate v 
from the equations of n and c. 

The hydrodynamic equations for active filaments (un- 
dergoing contractility and treadmilling) can be written 
as jH], 

<9 t n + A(n-V)n = K x V 2 n + X 2 V(V ■ n) + CVc 

+an- / 3|n| 2 n + f (1) 
<9 t e=-V-J = - V • Oo c n - DVc) (2) 

to lowest order in gradients and fields (contributions from 
v appear at higher order). The right hand side of 
and |2| represent contributions to active forces/torques 
and current J. The parameter a measures the deviation 
of the mean filament concentration c from the Onsager 
value which fixes the transition to orientational order. 
We choose a value of a and /3, such that the magni- 
tude of n in the ordered phase is close to unity in most 
places. The terms vq and A are uniquely active in origin, 
and represent an active advection and a nonlinear active 
convective contribution, respectively |12| . The athermal 
noise f is taken to be white with zero mean and variance 
equal to Ta/c, where Ta is the active temperature. We 
have dropped the additive noise term in the c equation, 
since the multiplicative noncquilibrium driving from the 
first term in J is more dominant. 

The values of these parameters depend on the micro- 
scopic active processes controlling actin dynamics. In 
a typical cellular context there are several microscopic 
active processes occurring simultaneously, such as tread- 
milling and actomyosin contractility. We therefore inde- 
pendently vary the parameters entering the dynamical 
equations over a range of values. Note that for the con- 
tractile motor-filament system, £ > 0, which is opposite 
to the flocking case considered in [TT| . With these param- 
eters, one can construct the following independent length 
scales - (i) correlation lengths L c , given by yjKi/a and 
\J K2 1 a, (ii) extrapolation lengths L e , given by Ki/Q and 
K2/C an d (hi) Peclet length, L p = D/vq, the ratio of dif- 
fusion coefficient to advection. We will predominantly 
work in the regime where both L c and L e are small, fur- 
ther, we study the phase diagram when the magnitude 
of A is zero or small. 

We first explore the phase diagram when the active 
temperature, Ta — 0. We find that the uniform orien- 
tationally disordered phase is stable when a < and 
C > but small, but beyond a threshold ( = —Da/vo, 
there is a clumping instability whose scale is set by the 
inverse of the maximal unstable mode kd — ^5/K\D, 
where 6 — C^o + Dot. On the other hand, the uniform 
orient ationally ordered phase which is stable when a > 
and £ = 0, is spontaneously unstable to splay distortions 
as soon as £vo > 0. Taking the ordering direction to be 
along x, we find that there is a band of unstable wavevec- 
tors centered around k x — 0, k y = k = (d+Ki+K z)' 

The final steady state configurations when a and £ are 
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Figure 2: (a) Logarithmic plot (base e) of the peak ampli- 
tude of the aster-density structure factor S evaluated at the 
reciprocal lattice vector G of the square crystal, as a func- 
tion of system size L, over a range of active temperatures 
T A (displayed symbols to the right). For < T A < 160, 
this amplitude scales as L 2 , suggesting true long range posi- 
tional order. Beyond this temperature, the amplitude scales 
as L 2 ~ v , where the critical exponent 77 is a function of Ta (in- 
set), suggesting a transition to a phase with QLRO. At still 
higher active temperatures, there is a transition to a short 
range order (SRO). (b) At T A = 240, corresponding to the 
QLRO phase, the structure factor evaluated about the peak 
{q\\ is parallel to G) exhibits a finite-size scaling form. The 
different system sizes L are shown as symbols. Inset shows a 
logarithmic plot (base e) of the scaling of the peak amplitude 
of 5* with L with a value of r\. Here, £ = 100 and the rest of 
the parameters as in Fig. 1. 

positive, depend on the extrapolation length L e relative 
to the other lengths, and can be obtained by numeri- 
cally solving Eqs. ([!]), Q. For this purpose, it is conve- 
nient to convert the equations to dimensionless form, by 
choosing the units of length, time and |n| to be D/2.5vo, 
D/2.5vq and respectively. The values of the var- 

ious parameters chosen for the numerics are written in 
these units. We use an implicit alternate direction opera- 
tor splitting scheme with length and time discretizations 
chosen to be Ax = 1 and At = 0.01 respectively |13| . 
with initial conditions for c(r,t) and n(r,t) being ho- 
mogeneous and random. The boundary conditions are 
chosen to be periodic; however, our results hold for other 
boundary conditions as well, as long as the system size 
L L c (when boundary effects negligible). In order 
to maintain conservation and non-negativity of the lo- 
cal concentration, we use symmetric spatial derivatives 
(that add up to zero over the whole system) and choose 
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Figure 3: (a) £-Ta phase diagram where the phases are (i) 
Boojum stripes (ii) Aster street (iii) Aster Lattice (iv) Bond 
tetratic (v) Bond Nematic and (vi) Isotropic phases, (b) Vari- 
ation of the bond orientational order parameters |^2| (ne- 
matic, brown) and |^4| (tetratic, blue) with Ta across the 
(iv) — > (v) transition. Parameter values as in Fig. 2. 



an adaptive (small) grid size At. 

Since ( > 0, the steady state configurations generically 
consist of a collection of defects such as Boojums, inward- 
pointing asters |14[ 115] . inward-pointing spirals or walls. 
This might be expected, since the filament current J oc n 
and the steady state equations for n is roughly a vector 
Poisson equation with ^j-Vc as source. 

We have explored the small extrapolation length 
regime in some detail, where upon increasing £, we en- 
counter the following defect phases (Fig. 1) : (i) Boo- 
jum stripes - consisting of alternate stripes of filaments 
oriented along x and configurations associated with a 
+2 defect known as "boojum" configurations with scale 
1/fco- (ii) Aster street - consisting of alternate stripes 
of filaments oriented along x and inward-pointing asters 
with scale 1 /ko and an aster size given by L p , the Peclet 
length, and (iii) Aster lattice - consisting of a square 
lattice of inward-pointing asters where the aster size is 
again L p and the 'lattice spacing' is l/k (> L p ). 

To see why asters should settle into a square lattice, 
note that in the limit A = 0, the right side of (IJ can 
be written as a derivative of an 'energy-functional', thus 
steady state solutions are minimisers of this 'energy' - 
we find that the 'energy density' of asters (in units of 
K\) arranged in a square (E sq = —35) is lower than in 
a triangular (E tr = —29) unit cell. As a check, we have 
verified that this square lattice persists when we reduce 
the spatial discretization to be much smaller than the 
aster lattice spacing. 

Since (n) ^ implies a movement of the active fila- 
ments with respect to the medium, the Boojum stripe 
and aster street are moving phases, while the aster lat- 



tice is a stationary phase. The mean drift velocity given 
by | (cn) | shows a discontinuous jump across the Boojum- 
aster street transition and a continuous transition at the 
street-lattice phase boundary (Fig 1(B)), where we find 
|(cn)| - |C - Cc| 7 with 7 w 1.38 ± 0.05. While we do not 
have an analytic explanation for the exponent value, it is 
easy to see why 7 > 1 - as one approaches the aster lat- 
tice phase from the street side, the filaments are drawn 
into the asters from the nearby parallel filaments by the 
£Vc term, as a result of which there is a nonlinear pos- 
itive feedback which draws in filaments more strongly, 
enhancing the rate at which the net drift velocity van- 
ishes. 

Other defects such as vortices or outward-pointing 
asters are unstable, since the current J oc cn ([2j. For 
instance, a vortex described by n cx eg is unstable to 
radial fluctuations with a rate ~ £vo/R, where R is the 
radius of the vortex. To stabilize such vortex configu- 
rations, one would need to include the binding to cross- 
linking proteins. One may obtain inward spiral asters, 
when the bend and splay distortions are comparable or 
even when the values of A are large. Being active such 
spiral asters would rotate with an angular velocity [5]. 
The sequence of defect configurations (Boojum — > aster 
— > spiral) obtained here is similar to the sequence of achi- 
ral tilt textures in a circular domain obtained on energy 
minimization |16l ITT] . In a different regime, when the 
extrapolation length L e is large, one obtains transient 
configurations of moving walls with filaments oriented 
normal to it. This resembles the configurations seen in 
reconstitution experiments [5J [3] and in recent numerical 
simulations [BJ. 

We now study the effects of active stochasticity on the 
steady state actin patterns by numerical integration of 
Eqs. ([TJ , ^2J) with noise, whose strength is parametrized 
by the active temperature, Ta > 0. At high £, where the 
steady state is an aster (square) lattice at Ta = 0, low 
noise results in phonon vibrations of the aster lattice. To 
quantify the state at low Ta, we compute the structure 
factor S(q) from the fluctuations of the density of asters, 
which we define as p(r) = — c(r)V • n(r). The structure 
factor shows Bragg peaks indexed by the reciprocal lat- 
tice vectors of a square, the amplitude of the peaks scales 
as L? for < T A < T* A (Fig. 2) - unlike the 2d equilib- 
rium solid, the 2d active solid shows true long range order 
(LRO) at finite T A ! Beyond T A , the amplitude of the 
Bragg peaks scales as L 2 ~ v ^ Ta \ suggesting a transition to 
a solid with quasi-long range order (QLRO). We confirm 
this from a scaling plot of the structure factor S(q) versus 
qL 2 ~ n for g|| G, the reciprocal lattice vector G = [1,1] 
(Fig. 2). In addition we compute bond-orientational or- 
der parameters, namely tetratic ^4 = (e l4e ) and nematic 
^2 = (e l2e ), where 9 is the orientation of the bonds be- 
tween nearest neighbour asters and the x-axis. As seen 
in Fig. 2, the bond orientational order parameters ^4 and 
^2 clearly indicates a discontinuous transition from a 
QLRO solid with tetratic order to a bond-nematic liq- 
uid. This jump in the order parameter $2 decreases on 
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Figure 4: (a) Sharp drop in density of filaments within asters 
(log plot in blue) with Ta and the corresponding variation 
of the order parameter cn (red). Parameter values same as 
Fig. 2. (b) Logarithmic plot (base 10) of the correlation func- 
tion C(q) = (n(q, t)-n(— q, t)) vs |g±|L for the ordered (solid) 
and exponential phases (circles) for different system sizes. 
The black solid line shows fit to power law C(q) ~ q~° with 
6 — 1.187 ± 0.061, which compares well with the Toner- Tu 
exponent = 6/5. Arrows show the lower (upper) wavevec- 
tor cutoffs due to system size (inter-aster distance). Values of 
£ = 10 and Ta = 100, rest same as Fig. 1. 

decreasing £ and approaches zero at a multicritical point. 

These structural transitions driven by the activity tem- 
perature Ta are associated with a remodeling of the 
asters; beyond the active solid phase the asters break and 
reappear transiently, resulting in a decrease in the mean 
aster density with increasing Ta [19] . The distribution of 
lifetimes of the asters is a power-law in the tetratic phase 
(-P(t) ~ t -2 7 ) and exponential in the nematic phase 

Interestingly, the decrease in the number density of 
asters is accompanied by an increase in the net polar- 
ization (cn) which jumps from being zero in the solid 
to a nonzero value in the tetratic and nematic phases 
(Fig. 2). Thus the effect of the active noise is to increase 
the polar order, weaning away filaments from the asters. 
Simultaneously, we find that orientational fluctuations 
about the ordered direction are massless. Figure 3 shows 
a power-law fit to the (unconnected) correlation function 
C(q) = (n(q, t) ■ n(— q, t)) versus q± (wave-vector per- 
pendicular to the ordering direction, taken to be along x; 
deviations occur both at low q± (corresponding to sys- 
tem size L) and high (corresponding to distance between 



asters, which decreases with increasing £). For compari- 
son, we plot the form expected when Q < 0; in this case, 
an exact RG calculation demonstrates that C(q) = q ± 6 ^ 5 
in d = 2 [IT]. The agreement with our numerics suggests 
that even for small, positive values of £, when the fila- 
ments are focussing, the orientational fluctuations in the 
high noise regime are controlled by the Toner- Tu fixed 
point. With a further increase in Ta, the system eventu- 
ally settles into an orientationally disordered phase. 

The sequence of transitions described above was in the 
high £ regime. At lower values of (, when the zero tem- 
perature phase is an aster street, an increase in Ta drives 
the steady state directly into a bond nematic phase. De- 
creasing £ further, i.e., starting from the Boojum stripe 
phase, we find that an increase in Ta leads to elliptical 
domains of Boojums whose scale is set by the correla- 
tion length, L c and D/vq. In fact, even when £ < and 
small, a nonzero Ta produces these elliptical domains of 
Boojum, as a consequence of the advection current vo. 

Within the calculation presented, the scale of the asters 
is set by the ratio D/vq, though simple extensions of 
Eqs. ([l|,([2j to include higher order terms or other pro- 
cesses can lead to an increase in the aster scale. One 
simple way is to increase the local concentration of fila- 
ments so as to attain a jammed aster configuration due 
to steric hinderance. Indeed the same result can be ob- 
tained with a lower concentration of filaments with the 
help of crosslinkers. An interesting alternative is to have 
a higher depolymerization rate kd at the aster cores bal- 
anced by a polymerization at the periphery. This would 
make the core size and hence the aster size bigger. A 
dramatic example of this phenomena is in the immuno- 
logical synapse, where large actin asters, of order several 
microns, drives the clustering of T-cell receptors [20J. 

To conclude, we have made a detailed study of the 
nonequilibrium phases and transitions of active filaments 
in two dimensions, both with and without active noise. 
The collective behaviour of active filaments gives rise to a 
variety of defect phases, which should be experimentally 
observable. We thank S. Mayor, S. Ghosh, P. Srivastava, 
S. Ramaswamy and S. Sengupta for discussions and a 
critical reading of the manuscript. MR acknowledges re- 
search grants from HFSP and CEFIPRA-35104. 
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